Spin generation away from boundaries by nonlinear transport 
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In several situations of interest, spin polarization may be generated far from the boundaries of 
a sample by nonlinear effects of an electric current, even when such a generation is forbidden by 
symmetry in the linear regime. We present an analytically solvable model where spin accumulation 
results from a combination of current gradients, nonlinearity, and cubic anisotropy. Further, we show 
that even with isotropic conductivity, nonlinear effects in a low symmetry geometry can generate 
spin polarization far away from boundaries. Finally, we find that drift from the boundaries results 
in spin polarization patterns that dominate in recent experiments on GaAs by Sih et al. [Phys. Rev. 
Lett. 97, 096605 (2006)]. 



Spin polarization can be generated and manipulated 
in semiconductors by means of electric fields and spin- 
orbit coupling. A prominent example is the spin Hall 
effect LLEl 3 l 4 l 5 l 6 l 7 l 8 El where a homogeneous electric current 
passing through a sample induces spin polarization s(r) 
near lateral edges, with opposite polarization at opposite 
edges. For rectangular homogeneous samples, this spin 
polarization falls off exponentially when moving away 
from the boundary on the length scale of the spin dif- 
fusion length L s , so for samples of size L > L s no spin 
polarization due to the spin Hall effect is expected far 
away from the edges (on scale L). However, in recent 
experiments on low-symmetry samples, Sih et aZ.^ ob- 
served polarized spins away from the edges of a GaAs 
sample subjected to an electric current and concluded 
that there are transport effects beyond the simplest spin 
Hall effect near edges. 

In existing analytical theories of the spin Hall effect, 
spin polarization in spin-orbit coupled media has been 
considered in the linear transport regime for small electric 
field E. Here we show that in this regime, spin generation 
away from boundaries is forbidden or strongly suppressed 
for sufficiently symmetric samples. However, in the con- 
text of extrinsic effects, by considering nonlinearities in 
charge transport , we find a new mechanism of generating 
electron spins by electric current. This nonlinear regime 
is of practical importance, because experiments are of- 
ten performed in a range of elect ric fi elds with nonlin- 
ear current-voltage characteristics We then present 
an analytically solvable model in which the bulk spin 
generation in a radially symmetric sample geometry is 
due the nonlinearity and the anisotropy of conductivity 
tensor. Finally, we numerically solve the charge trans- 
port and spin drift-diffusion equations for the sample ge- 
ometries used in Ref. [TOj The patterns of spin accumula- 
tion we find strongly resemble experimental findings. We 
establish the existence of two contributions to the spin 
polarization: the first contribution is generated at the 
boundaries and then drifts large distances (while diffus- 
ing further away from the boundaries), and the second 
one is generated away from boundaries. While we find 
that in the recent experiments the first of these domi- 
nated, we propose setups that should allow unambiguous 
observation of the spins generated in the bulk. 



We consider a diffusive system with a weak extrinsic 
spin-orbit coupling and typical system size L. The spin 
current contains drift and diffusion contributions. Ad- 
ditionally, electrical current, via the extrinsic spin-orbit 
coupling, induces a spin Hall current j SH , which acts as 
a source for the spin polarization s z (r) with a genera- 
tion rate T z = — divjg H . For length scales much larger 
than the mean free path £, spin-orbit processes (e.g., 
Dyakonov-Perel or Elliot -Yafet mechanisms) lead to a fi- 
nite spin lifetime t s . Together with a spin diffusion coef- 
ficient D s , the spin lifetime defines a spin diffusion length 
L s = V D s t s . To analyze the spin polarization away from 
boundaries, we consider the regime £ <C L s < L. On 
length scales large compared to I and in the absence of 
an external magnetic field, the spin density s z (r) obeys 
the drift-diffusion equation 

s z = div (D s Ws z ) - div (v dr s z ) + T Z (1) 

where drift velocity Vd r is proportional to the local elec- 
tric current density. This introduces yet another length 
scale, the drift length — Vd r T s . 

In the linear transport regime, one can evaluate T z 
by writing the spin current induced by the electric field 
as j| H = <r SH E, with the spin Hall conductivity ten- 
sor <x that does not depend on E. Assuming nonin- 
teracting electrons, in the absence of the intrinsic spin- 
orbit coupling, one can independently consider the spin 
species with opposite polarization and then relate the 
spin Hall effect to the anomalous Hall effect.' In par- 
ticular, at magnetic field B — 0, we find da /dep = 
(h/2^e)da Aii /dB\ B =o, with o- AH = \ [a(B) - a(-B)], 
and where a(B) is magnetic field-dependent charge con- 
ductivity, £p is the Fermi energy, fi is the magnetic mo- 
ment, and e is the electron charge. From the Onsager re- 
lation for a(B), we see that <r SH must be anti-symmetric. 
Using this property for a homogeneous two- or three- 
dimensional system, we find 

r - = ~\ E« fc 4 He «fe(V x E) fc = 0, (2) 

i.e., the spin generation rate vanishes and no spin polar- 
ization would be found far away from boundaries in the 
linear regime. 
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For a system with intrinsic spin-orbit interaction, the 
interpretation of th e spin currents defined above be- 
comes less clear 12HS] Thus, instead we consider the elec- 
trically induced spin density s(r) inside the sample on 
length scales larger than L s . For typical sample dimen- 
sion L > L s , the field E changes slowly over L s , and 
we make a gradient expansion linear in E, leading to 
Si(r) « J2j AjEj(r) + J2 jk B ijk djE k {r). For GaAs, the 
cubic symmetry group contains mirror planes, im- 
plying that s transforms as a pseudo- vector; thus all 
Aij vanish. The only pseudo-vector linear in V and E 
is V x E = 0. Going to the next order in the gradi- 
ent expansion, there are contributions cx (9~ x — d 2 y )E z 
and oc d z {d x E x — d y E y ), where x, y, z are the principal 
crystallographic axes. These terms also vanish in the in- 
plane geometry (d z Ei = 0) with an in-plane electric field 
(E z = 0). 

We have found that spin generation away from bound- 
aries may be forbidden or strongly suppressed in the lin- 
ear transport regime. Therefore, we will now discuss how 
nonlinear effects can lead to such spin generation. We will 
restrict ourselves to the case of extrinsic spin currents 
that are generated by spin-dependent impurity scatter- 
ing; they are given by j| H = z x (^J c — 2An|E)P where 
the two terms correspond to the skew-scattering and side- 
jump contributions, respectively. Here, J c — ctE is the 
charge current, but now we allow the conductivity ten- 
sor to be a function of the electric field. Further, A is a 
material-dependent spin-orbit coupling constant and the 
skewness 7 is of order A but it also depends on the prop- 
erties of scatterers and on electron distribution function. 
The side-jump part of the current does not contribute to 
T Zl and we obtain 
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Assuming that 7 is constant, that the conductivity is 
isotropic, and that its position dependence is described 
by the local electric field, cr(r) — a(E(r)), Eq. ([3| sim- 
plifies to 



(4) 



In samples where the extrinsic spin Hall effect was 
observed 3 -^ th e charge conducti vity a(E) was found to 
be an increasing function of fieldpES we expect that T z 
is finite in these samples, and spin polarization can be 
generated away from boundaries (sec discussion below) . 

Next, we show how an anisotropic, nonlinear conduc- 
tivity leads to spin generation away from the boundaries. 
In general, according to Eq. pi, inhomogeneous electric 
fields are required for a finite T z . Thus, we now analyze 
the inhomogeneous field in a Corbino geometry, where 
a total current I is injected at r — a into an infinite 
two-dimensional sample. In a (001) film of crystal of full 
cubic symmetry, the leading nonlincarities in J C (E) are 



where the components i = x,y are taken along the prin- 
cipal crystal axes and we assume that ai^E 2 (a) <C a. 
One can expect anisotropic nonlinear terms of consider- 
able magnitude for many-valley semiconductors like SiGe 
quantum wells, and also for Al^Gai-^As quantum wells 
with x close to the direct-indirect gap transition. We first 
solve for the electrostatic potential 3>(r); in polar coor- 
dinates it will be of the form ^ m <I> m (r) cos 4mcp. The 
term $o( r ) will not contribute to T z [Eq. pi], so we 
next consider the lowest harmonic, 3>i(r) cos4i/j. Because 
div J r = 0, we see that 
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Requiring that $i(r — > 00) = 0, we obtain 
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The spin generation rate is then 
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Jl = (a + a 2 E 2 )E l + a 1 E?. 
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So indeed, the combined anisotropy and nonlinearity of 
conductivity lead to a spin generation, which, for r 3> 
L s , results in the spin density is s z (r, tp) = T z (r, ip) t s cx 
7 3 sin(4c/3) /r 4 . We emphasize that this polarization falls 
off only as a power law. Furthermore, it consists of four 
sectors of up-spins separated by four sectors of down- 
spins. 

We now analyze the spin polarization s z (r) in sys- 
tems with isotropic (nonlinear) conductivity but with 
less symmetric geometries. We consider experiments on 
unstrained GaAs samples by Sih et al.. 10 The spin Hall ef- 
fects in s uch s amples is believed to be primarily extrinsic 
in origin G3M1 j n these experiments a T-shaped geometry, 
as shown in Fig.[l] was used. In an electric field, electrons 
flow from the bottom to the top of the main channel, 
some of them entering into the side-arm; the gradient 
of the electrical field E(r) becomes large in the region 
near the entrance of the arm, and a field-dependent con- 
ductivity will then lead to spin generation across this 
region. Furthermore, due to the spin Hall effect, spins 
are also generated near the sample edges and can then 
diffuse and drift along the electric field into the center of 
the side arm — below, we find that this latter mechanism 
dominates the experimental observations. We choose re- 
alistic values of the parameters as follows. Electron den- 
sity n = 3 x 10 16 cm -3 and sample dimensions are taken 
from Ref. 1101 Further, we assume that the sample is ho- 
mogeneous and the spatial dependence of a, r s , and D s is 
controlled by the local field E(r). Low field conductivity 
c(E) for GaAs samples was obtained from unpublished 
datap^ and for higher fields we took the £ , -dependcncc 
of the conductivity measured for Ino.07Gao.93As in Ref. 3 
as a guideline, which increased by a factor of two when 



the field increased from to 20 mV//jm. The spin relax- 
ation time t s (E) was loosely based on the experimental 
data taken at two points away from the boundaries for a 
range of electric fields^ we show our assumed o{E) and 
t s (E) in Fig. 1(c). Furthermore, it was found in Ref. 3 
that the spin diffusion length L s was field-independent 
within error bars; thus we take D S (E) — L 2 s /t s (E) with 
constant L s . We chose L s = 7 fim found from the best 
fit to the data of Ref. [TH] for the main channel. Finally, 
we take spin-orbit coupling co nstan t A = 5.3 A 2 and es- 
timate the skewness 7 = 1/700.^™ In our simulation, we 
first solve for the electrostatic potential ^(r) for an ap- 
plied dc voltage and determine the spin generation rate 
T z (r) [Eq. ^] and the drift velocity Vd r (r) = J c (r)/ne. 
Because there are no direct indications of a considerable 
spin relaxation at the boundaries, we consider the spin- 
conservation boundary condition n ■ (D s Ws z — jg H ) = 0. 
This accounts for the spin generated at the boundaries 
due to the spin Hall effect. We then solve the spin drift- 
diffusion equation [Eq. and find the stationary spin 
polarization s z (r). 



The simulated s z (r) is shown in Fig.JTJa), taking exper- 
imental values for sample dimensions and electric fieldP^ 
The most striking feature is a spin distribution inside the 
side-arm: it results from spins generated at the bound- 
aries near the corner of the lower part of the side-arm, 
which diffuse and then drift a distance Ld r = Vd r T s , which 
is much longer than L s . Conversely, near the upper part 
of the side-arm the spin population is very low, since 
spins generated at the boundary drift out of the side- 
arm; i.e., the strong drift (Ldr > L s ) leads to an asym- 
metric spin distribution. However, in Ref. [10] a square 
wave voltage V was applied for lock-in detection; then 
the measured spin polarization, | [sz(V) — s z (—V)], be- 
comes symmetric, and we show the 'symmetrized' po- 
larization in Fig. |TJb) . We find good agreement with 
experimental data, see Fig. 1(c) in Ref. [TO] in particular, 
the maximum of spin distribution is at a similar position 
inside the side-arm, about 10 - 15 fim deep. In Fig.[TJc), 
we show s z along the central section of the side-arm for 
different depths of the side-arm and find that for depths 
of 10 or 20 //m, the spin population is maximal at the 
right edge of the arm. 

The simulation includes contributions from spins gen- 
erated away from boundaries with rate T z [Eq. Q3p] . For 
comparison, we show the spin polarization for T z = 0, 
see dashed line in Fig. 1(c). We, however, find that the 
absence of T z changes the resulting spin polarization by 
less than 15%. Therefore, in the parameter regime of 
the experiments of Ref. 1101 the dominant contribution 
to the stationary spin distribution inside the side-arm 
comes from the drift of the spins generated at the sam- 
ple boundaries. This results from the large drift distance 
Ld n which is of the same length scale as the width of the 
side-arm. 
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Figure 1: Simulated spin accumulation s z (r) for experimen- 
tal geometry and parameters of Ref. 1101 Electron flow is 
from bottom to top in the main channel, and spills over into 
the side-arm. Distances are measured in microns and elec- 
tric field in the main channel is 9.5 mV/(im. (a) Numeri- 
cal solution of spin drift-diffusion equation [Eq. |TJ] for a 40 
/im sidearm in the dc regime, (b) Symmetrized spin accu- 
mulation | [sz(V) — s z (— V)]. (c) Spin polarization s z (x,y) 
(convoluted with a Gaussian with standard deviation of 1/im 
that is associated with laser spot) for y = (along dashed 
line in (b)) and for side-arms of depth 10 (green), 20 (red), 
and 40 /im (blue) is in good agreement with the experimental 
data in Fig. lc of Ref. [TO] The dashed line shows a simulation 
without bulk spin generation (Y z = 0) . The inset shows spin 
relaxation time r s (E) and charge conductivity cr(E) as used 
in our simulation. Panels (d), (e) show s z (r) on a different 
color scale for a sample larger by a factor a — 10. In (e), the 
spin generation rate Y z is neglected. 



In order to distinguish spins generated in the bulk from 
those originating at the boundaries, one needs a set up 
where both L s and L^ r are smaller than the width of the 
side arm, L. While L s < L was satisfied in the experi- 
ments of Ref. [TOland in the simulations above, < L 
was not. The inequality can be achieved (i) by reducing 
Ld r with a smaller applied field (and, therefore, smaller 
drift velocity), (ii) by reducing Ld r by pulsing the ap- 
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plied fields and thus allowing drift only for a shorter time, 
(iii) by using an alternating sequence of pulsed electric 
fields, which can largely cancel the drift and the spin gen- 
eration at the boundaries, or (iv) by increasing sample 
dimensions L. This exponentially suppresses the contri- 
butions from the boundaries (as long as L s is sufficiently 
short), making the relative contribution of the genera- 
tion away from boundaries arbitrarily large. However, 
this also suppresses the absolute value of spin polariza- 
tion, as we discuss now. 

For approach (i), a smaller drift velocity is achieved 
by choosing a sufficiently small electric field, so that 
Vdr <C L/t s . Because of the scaling T z oc a'(E)E 2 /L, 
with cr'(E) = da/dE, small drift velocities also im- 
ply that spin polarization due to generation away from 
boundaries is weak, (ii) When the electric field is only 
applied for pulses of duration t p , which are short com- 
pared to r s , the drift length is reduced to Ld r = i P Vd r - 
Also, the pulses must be at least r s apart, thus the aver- 
age spin polarization signal is reduced from the dc case 
by a factor smaller than i p /r s . For example, applying 
electric field pulses with t p = 1 ns and separation time 
20 ns, for the parameters used in Fig. [T] spin generation 
away from boundaries could be observed with a time- 
averaged value of s z < 0.2 /mi~ 3 . Another approach (iii) 
is to use alternating pulses of large positive E for dura- 
tion t\ and negative E of smaller magnitude for a longer 
time t2, such that the time-averaged electric current is 
zero. Then, if the period t\ + ti is smaller than both t s 
and £/vd r , the effects of spin drift should be largely can- 
celed. Furthermore, during t\ and equal amounts of 
spin polarization are generated at the boundary, but with 
opposite signs, which leads to a cancellation of the contri- 
bution of spins generated at the boundaries. In contrast, 
because the bulk spin generation T z is nonlinear, it can 
have a nonzero time average. Thus, the effectiveness of 
this scheme will depend on the strength of nonlinearitics 
involved. 

Note that in approaches (i)-(iii), diffusion of spins gen- 
erated at the boundaries may still dominate s z (r) in the 
bulk. For example, in simulations with the geometry of 
Fig. [T]but for smaller fields, the bulk s z still contains a 
large relative contribution from spins generated at the 
boundaries. To eliminate such a diffusion effect, one 
could use approach (ii) in a pump-probe scheme, where 
the spin polarization is detected shortly after the pulse, 
providing a direct measurement of T z (r) . 

The most straightforward way to reveal spin generation 
away from boundaries is to increase the sample's linear 
dimensions by a factor a big enough so that the sidearm 
opening exceeds the drift length Ld r , while keeping the 



average electric field in the main channel fixed (iv) . This 
reduces T z and hence the bulk spin polarization by a fac- 
tor a. However, the signal to noise should not decrease 
since the signal can be taken from an area that is larger by 
a factor of a 2 , thus reducing the noise by a. Figs. l(d,e) 
demonstrate spin populations in a sample 10 times big- 
ger than the sample of Fig. [T] A noteworthy feature is 
found in the vicinity of the corners: the spin population 
there is of the opposite sign than the population on the 
adjacent boundary. This is a distinctive feature of spin 
generation away from boundaries. For comparison, we 
show in Fig. 1(e) the spin polarization for T z = 0. No 
spin polarization is found in the central part of the side- 
arm; therefore, the spin polarization in the central part of 
Fig. 1(d) indeed results from spins generated away from 
boundaries. 

In conclusion, we have shown that, in a linear trans- 
port regime, the generation of spin polarization away 
from boundaries is forbidden or strongly suppressed. 
However, in nonlinear transport, spins can be generated 
away from boundaries and we analyze such generation 
resulting from spin Hall currents due to the "extrinsic" 
skew scattering mechanism. For an anisotropic nonlin- 
ear charge conductivity tensor, we analytically evaluate 
the spin generation in a Corbino geometry. We also sim- 
ulate the spin polarization in a T-shaped geometry us- 
ing isotropic nonlinear conductivity and field-dependent 
spin lifetimes, appropriate to recent experiments. We 
find that the spin accumulations in the experiments of 
Ref . [TO] were primarily due to drift of spins generated at 
the boundaries, but we suggest other setups, where drift 
should be relatively unimportant, and nonlinear spin gen- 
eration away from boundaries should be observable. 
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PHYO 1-1 7795, and by Harvard Center for Nanonscale 
Systems. 

Note Added: In a recent work, Stern et aZP^^ have 
independently shown that drift alone can account for ob- 
servations of Ref. Uni Golizadeh-Mojarad and DattaP^ 
have obtained spin polarizations qualitatively similar to 
those seen in experiments^ using a model with intrinsic 
(Rashba) spin orbit coupling. However, the spin coher- 
ence length and sample dimensions used were two or- 
ders of magnitude smaller than the experimental values. 
Pershin and Di VentraP^ considered rectangular samples 
with inhomogeneous charge density n, which leads to 
spin generation linear in the electric field, assuming that 
a cx n [cf. Eq. (^J]. 
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